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An improved Hamiltonian constraint operator is introduced in loop quantum cos- 
mology. Quantum dynamics of the spatially flat, isotropic model with a massless 
scalar field is then studied in detail using analytical and numerical methods. The 
scalar field continues to serve as 'emergent time', the big bang is again replaced by 
a quantum bounce, and quantum evolution remains deterministic across the deep 
Planck regime. However, while with the Hamiltonian constraint used so far in loop 
quantum cosmology the quantum bounce can occur even at low matter densities, 
with the new Hamiltonian constraint it occurs only at a Planck-scale density. Thus, 
the new quantum dynamics retains the attractive features of current evolutions in 
loop quantum cosmology but, at the same time, cures their main weakness. 
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I. INTRODUCTION 

The spatially flat, isotropic model was recently investigated in detail in the setting of loop 
quantum cosmology (LQC) jH. (For a brief summary of results, see 0.) That investigation 
introduced a conceptual framework and analytical and numerical tools to construct the 
physical sector of the quantum theory. These methods enabled one to systematically explore 
the effects of quantum geometry both on the gravitational and matter sectors and to extend 
the previous results in LQC. 

If the only matter source is a massless scalar field, every classical solution has a singularity 
which, furthermore, persists in the Wheeler-DeWitt theory. Therefore, to bring out effects 
which can be unambiguously attributed to the quantum nature of geometry underlying 
LQC, this model was analyzed in detail in 0, The analysis led to the following results: 
i) the scalar field was shown to serve as an internal clock, thereby providing a detailed 
realization of the 'emergent time' idea; ii) the physical Hilbert space, Dirac observables and 
semi-classical states were constructed rigorously; and, iii) the Hamiltonian constraint was 
solved numerically to show that in the backward evolution of states which are semi-classical 
at late times, the big hang is replaced by a quantum bounce. Furthermore, thanks to the non- 
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perturbative, background independent methods, unlike in other approaches the quantum 
evolution is deterministic across the deep Planck regime. 

These results are attractive and add to the growing evidence B 0] suggesting that the 
quantum geometry effects of loop quantum gravity (LQG) hold a key to many long stand- 
ing questions. Control over the physical sector of the theory and availability of numerical 
methods also enabled to us to analyze the physics of the bounce in greater detail than could 
be done in previous investigations. These details enriched our understanding of both the 
semi-classical and deep Planck regimes. However, they also brought out a serious limitation 
of the framework: the critical value pcrit of the matter density at which the bounce occurs 
can be made arbitrarily small by increasing the momentum of the scalar field (which 
is a constant of motion). Now, large values of are permissible and even preferred for 
semi-classical considerations at late times. Since it is physically unreasonable to expect 
quantum corrections to significantly modify classical predictions at low matter densities, 
this limitation is a serious drawback. Indications of such problems had appeared also in 
some of the earlier results based on effective equations (see, e.g., However, since the 

approximations inherent to effective equations break down in the deep Planck regime, one 
did not have an a priori reason to believe that these were consequences of the underlying 
quantum equations of LQC, rather than artifacts of approximation schemes. The analysis 
of in unambiguously showed that the problem lies with the quantum equations themselves. 
Thus, while the quantum evolution predicted by the existing LQC constraint operator has 
a number of attractive features, it has one major weakness. A key question then is whether 
one can change the definition of the constraint operator in a subtle way so that this problem 
is overcome, but the attractive features of the evolution are all preserved. 

As reported in P, 0], the answer is in the affirmative. The purpose of this paper is 
to justify this assertion by providing the detailed construction of the new Hamiltonian 
constraint and analyzing the resulting quantum dynamics. We will find that the required 
change is in fact physically well motivated and conceptually quite compeUing. 

The main ideas can be summarized as follows. Recall that, to obtain the expression 
of the quantum constraint, one has to first introduce an operator representing curvature 
of the gravitational connection. Now, a key feature of LQC — inherited directly from full 
loop quantum gravity (LQG) — is that while there are well defined quantum analogs of 
holonomies, there is no operator corresponding to the connection itself 0,11,0]. Thus, one is 
naturally led to define the curvature operator in terms of holonomies. In the classical theory, 
curvature can be expressed as a limit of the holonomies around a loop as the area enclosed by 
the loop shrinks to zero. In quantum geometry, however, one can not continuously shrink the 
loop to zero area; there is a smallest non-zero area eigenvalue, or an 'area gap 
The physical idea then is to incorporate the existence of the mass gap A 4] in the limiting 
procedure. However, in the existing definition of the quantum constraint, this physical idea 
was mathematically implemented in a rather indirect manner. Initially, it was meant to be 
only a first step to gain a qualitative understanding of dynamics within LQC. Nonetheless, as 
more and more physically appealing features of this dynamics were discovered, the constraint 
operator obtained by this procedure was taken more and more seriously and constituted the 
basis of a number of approximation schemes and effective descriptions. 

The new Hamiltonian constraint introduced in this paper is based on the same princi- 
ples but the curvature operator is now constructed by a more direct implementation of the 
underlying physical ideas. The structure of the final constraint operator does change sig- 
nificantly but, when suitably recast, its form is actually somewhat simpler. Furthermore, it 
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has additional conceptually attractive features. For example, its Wheeler-DeWitt (WDW) 
limit naturally comes with the 'covariant factor ordering' of that theory. Once the new 
constraint is introduced, we apply analytical and numerical methods of QJ] to extract phys- 
ical predictions from the theory. Specifically, we will build the physical Hilbert space Hphy 
from solutions to the new quantum constraint, introduce on it a complete family of Dirac 
observables, construct suitable families of states which are semi-classical at 'late times' (e.g., 
'now'), and numerically evolve them back in time. As in we will find that in this back- 
ward evolution, the big bang is again replaced by the quantum bounce. Qualitative features 
of the evolution are very similar to those observed in except that the big bounce occurs 
precisely when the matter density enters the Planck regime, irrespective of the value of Pif, 
and other initial conditions. 

To summarize, the Hamiltonian constraint will be improved by implementing the physical 
idea of j3| in a more satisfactory manner. The resulting dynamics will retain the attractive 
features of the older quantum evolution but free it from its main drawback. This rather 
delicate interplay between physics and mathematics is both interesting and instructive. We 
wish to emphasize however that the present procedure still retains a basic limitation of the 
existing treatments: the Hamiltonian constraint is not systematically derived from the full 
theory. Indeed, this important task can not yet be undertaken because there is no unam- 
biguous Hamiltonian constraint in the full theory which can serve as the natural starting 
point for a systematic reduction. 

The organization of this paper parallels that of and we use the same notation. We 
direct the reader to that paper for motivations as well as further technical details. In Section 
im we introduce the new Hamiltonian constraint operator. In Section IIIII we discuss the 
WDW theory that results when one ignores the effects of quantum geometry. As in P|, 
we find that the singularity is not resolved in that limit. In Section HVl we return to LQC 
and construct the physical sector of the theory. Quantum states which are semi-classical at 
'late times' are then numerically evolved backwards in Section We find that the classical 
big bang is replaced by a quantum bounce which occurs when the matter is compressed 
enough to acquire a density of the Planck scale. Thus, in the deep Planck regime, quantum 
geometry has the effect of making gravity strongly repulsive. A key virtue of the new 
constraint operator is that these effects are completely negligible when the matter density is 
significantly smaller than the Planck density. Section IVTl summarizes the main results and 
briefly discusses possible extensions. While our detailed analysis is restricted to a specific and 
rather simple model, the key ideas behind the construction of the new constraint operator 
can be applied much more generally. This is illustrated in Appendix ^ where we allow 
for the presence of a cosmological constant. (Other models will be discussed also by other 
authors in forthcoming papers.) Appendix El discusses conceptual issues that are especially 
relevant to approximation schemes used in constructing effective equations. 

A brief overview of LQG emphasizing these developments in LQC can be found in [13^ . 



II. THE IMPROVED CONSTRAINT OPERATOR 



This section is divided into three parts. In the first two we introduce the key new elements 
and in the third we put them together to construct the new Hamiltonian constraint operator. 
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A. Strategy 

Let us first collect the required background material and notation. (For details, see 
In a systematic Hamiltonian treatment of spatially fiat, isotropic models, one has to first 
introduce an elementary cell V and restricts all integrations to this cell. Fix a fiducial flat 
metric °qab and denote by Vo the volume of the elementary cell V in this geometry. The 
gravitational phase space variables — the connections A\ and the density weighted triads 
E'^ — can be expressed as: 

A\ = cV;''°uol, and Ef=pK~^v^°e^ (2.1) 

where (li;* , °e°:) are a set of orthonormal co-triads and triads compatible with °qab and adapted 
to the edges of the elementary cell V. Thus, the symmetry reduced gravitational phase 
space is only 2-dimensional, coordinatized by the pair (c, p). Thanks to the availability of 
the fiducial cell V, the coordinates (c, p) are insensitive to the choice of the fiducial metric, 
i.e., remain unchanged under the rescaling °qab — > a^°5'afe- The fundamental Poisson bracket 
is given by: 

{c, p} = (2.2) 

where 7 is the Barbero-Immirzi parameter. The gravitational part of the Hamiltonian 
constraint can be written as: 

Cg,av = -1'' [ d'xe^.ke-'E'^'E'^ = -^c"^ (2.3) 
Jv 7 

where e := A/|deti?|. 

In quantum theory, following Dirac, one first constructs a kinematical description. The 
Hilbert space H^j^ is the space L^(]RBohr5 c^/^Bohr) of square integrable functions on the Bohr 
compactification of the real line. To specify states concretely, it is convenient to work with 
the representation in which the operator p is diagonal. Eigenstates of p are labelled by a 
real number fi and satisfy the orthonormality relation: 

(/^i |/i2) = 5mi,/^2 • (2-4) 

Since the right side is the Kronecker delta rather than the Dirac delta distribution, a typical 
state in Ti.^^^ can be expressed as a countable sum; |\E') = J2n c*-"^|/^n) where c^"'^ are complex 
coefficients and the inner product is given by 

(^i|^2) = $^c^)4"). (2.5) 

n 

The fundamental operators are p and expiA(c/2): 

^ , , Sn^yip, . , , iXc , , , , , 

p|/i) = — - — and exp — l/i) = l/i + A) (2.6) 
o / 

where A is any real number. Since the holonomy h\^^ of the gravitational connection 
along a line segment A is given by: ^ 

= cosf I + 2sin^rfc (2.7) 



^ Here I is the unit 2x2 matrix and is a basis in the Lie algebra su(2) satisfying tj = ^EijkT^' — j^ij- 
Thus, 2iTk — Ok, where Ui are the Pauh matrices. 
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the corresponding holonomy operator has the action: 

hi'^ l/i) = ^ (l/i + A) + l/i - A)) I + ^ (l/i + A) - l/i - A)) Tfc . (2.8) 
However, just as there is no operator corresponding to the connection itself in full LQG 



riowever, mst as tnere is no operator corresponai 
0, i, i, ES, there is no operator c on 0, Q. 



/grav 



To describe quantum dynamics, we have to first introduce a well-defined operator on H^^j^ 
representing the Hamiltonian constraint Cgrav Since there is no operator corresponding to 
c itself, we return to the integral expression of the constraint from the full LQG: 

Cg,av = -7"' / d'xe,,k e-'E'^^E"^ F^, . (2.9) 



V 



For passage to quantum theory, we first need to express this classical constraint in terms of 
the elementar y v ariables p and h^,^^ which have unambiguous quantum analogs. As in the 
full theory |9[1l5l|. the term involving triads can be written as 

e-' E^^E'' = V o^abc o^k UW |/,(A)-i^ V} r,) (2.10) 

k 2tt^GXVo' ^ ' 

where V = is the volume of elementary cell V in the physical metric determined by 

p and the holonomy is evaluated along the segment A°e^ (i.e., a segment of oriented 
length A along the kth edge of the elementary cell V) . Note that this identity holds for any 
choice A, even when it is allowed to be a function of p. Let us allow for this possibility and 
fix the appropriate A at the end. 

For the field strength F^^, we use the standard strategy employed in gauge theories. 
Consider a square in the i-j plane spanned by a face of the elementary cell, each of whose 

1/3 

sides has length XVo with respect to the fiducial metric °qab- Then, 'the ab component' of 
the curvature is given by: 

F^, = -2 lim Tr " r'^ K • (2-11) 

Ara^O y X^Vo^^ Jab 

Here Ar\j is the area of the square under consideration, and the holonomy around the 
square Dij is just the product of holonomies (j2.7|) along the four edges of Dij: 

hS!^ = hi'^hf\h^Y\h^r. (2.12) 

Combining Eqs. ()2.10|) and 1)2.111) . Cgrav can be written as 
Cgrav = ^lim^C^^i^ where 



sin Ac 



ijk 
4 



87rG73 A 



sin Ac (2.13) 
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where in the last step we have used a symmetric ordering of the three terms for later con- 
venience. Since the constraint is now expressed in terms of elementary variables which have 
unambiguous quantum analogs, it is straightforward to write down the quantum operator 
Cgrlv. However, the limit Ar\j — > of this operator does not exist. This is not acciden- 
tal; had the limit existed, there would be a well-defined local operator corresponding to 
the curvature F*^ and we know that even in full LQG, while holonomy operators are well 
defined, there are no local operators representing connections and curvatures. This feature 
is intimately intertwined with the quantum nature of Riemannian geometry of LQG. The 
viewpoint in LQG is that the failure of the limit to exist is a reminder that there is an 
underlying quantum geometry where eigenvalues of the area operator are discrete, whence 
there is a smallest non-zero eigenvalue. A, i.e., an area gap. Thus, quantum geometry is 
telling us that, on physical grounds, it is inappropriate to let Ara go to zero; the 'correct' 
procedure must take into account the existence of the area gap. 

Up to this stage, we have followed the standard LQG reasoning, first introduced in j^. 
The difference comes in the implementation of the above idea. The viewpoint we now adopt 
is that since quantization of area refers to physical geometries, we should shrink the loop 
Dij till the area enclosed by it, as measured by the physical metric qab, reaches the value 
A. Since the physical area of faces of the elementary cell is \p\ and since each side of dij is 
A times the edge of the elementary cell, we are led to choose for A a specific function fi{p), 
given by 

fi^ bl = A = (2v^7r7)4i (2.14) 

While ft is a non-trivial function on the phase space, its analog fio in the existing LQG 
treatments was required to be a constant.^ Technically, this contrast turns out to be im- 
portant. If we write the quantum constraint using the \fi) basis, the constraint operator of 
[1, 4] is a difference operator with uniform step size, given by 4/Xo. In this basis the action of 
the new quantum constraint would be again given by a difference operator, but the step size 
now depends on the state it operates on. This difference turns out to be subtle enough 
to remove the major weakness of the current LQG constraint operator, while retaining its 
physically desirable features. Gonceptually, the new strategy has important ramifications 
for the fundamental curvature operator. In both treatments, it is non-local. However, while 
the non-local operator used in the existing literature depends only on the connection, the 
non-local operator introduced here depends both on the connection and the geometry. In 
the classical limit, effects of quantum geometry are negligible and the classical limits of both 
operators yield the standard expression of the curvature Fab. 

Finally, some heuristics involving full LQG can be used rather effectively to bring out 
the physical difference between the two strategies. In the sector of full LQG consisting 
of a small but a finite neighborhood of homogeneous isotropic cosmologies, one should be 
able to carry out partial gauge fixing, enabling one to speak of 3-metrics (rather than their 
equivalence classes under diffeomorphisms) . Let us consider quantum geometry states based 
on graphs and for simplicity focus on the fixed fiducial cell. Then, one would expect that, 
as the scale factor a grows, the number N of vertices in the graph contained in the fiducial 



^ The procedure used in LQC (and specifically in [lHjl) so far can be summarized as follows. One notes 
that, regarded as a state in the connection representation, each holonomy ft.^^""* is an eigenstate of the 
area operator (associated with the face of the elementary cell orthogonal to the fcth direction) . One fixes 
fio by demanding that this eigenvalue be A and finds /io = Sa/S/S, a constant. 
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cell would grow as = K^a^ for some constant Ko- In the Hamiltonian constraint of full 
LQG, the operator Fah would have to be evaluated at these vertices. A natural procedure 
is to introduce an 'elementary cube' around each vertex and calculate the holonomy around 
the faces of that cube. Since the number of elementary cubes contained in the fiducial 
cell increases as a grows, the area of their faces, as measured by the fixed fiducial metric 
would decrease as Ar°^^ ~ (1/A^)^/^ ~ Hence, as measured by the fiducial metric, 

the length of each edge of the face around which the holonomy is evaluated would go as 
1/ a/IpI — i.e. precisely as p,. Thus the functional dependence of /2 on /i can be thought of 
as the remnant left behind by the sector of the full theory of interest to cosmology, on the 
mini-superspace considered in this paper. Fixing the edge length to be /Xq, i.e., making it 
state independent, amounts to ignoring the 'creation of new vertices' that one expects to 
accompany the expansion of the universe. These considerations complement ideas recently 
advocated by Bojowald to incorporate inhomogeneities in LQC 



B. The operator expz(^c/2) 

Having fixed A to be the specific function fi{p), we are now led to write the expression of 
the quantum constraint by replacing A in the right side of ()2.13|) by ft, the holonomies and 
the volume functions by corresponding operators, and the Poisson bracket by 1/ih times 
the commutator. However, to carry out this task, we first need to define the operators 
cos(/ic/2) and sin(/ic/2), i.e. the operator expz(/ic/2). This task is a somewhat subtle 
because i) there is no operator corresponding to c and hence to flc, and, ii) since /i is a 
function of p, expi(/ic/2) can not be expressed as a function of the 'elementary variables' 
— p and expi{kc/2) where /c is a constant. 

However, in this task we are aided by geometric considerations. Let us begin with states 
^(//) in the Schrodinger Hilbert space, L^(M, d/i). For any real constant k we have: 

e'f ^(/i) = ^(/i + k) = e''^ ^(/i) . (2.15) 

Thus, the action of the operator is to drag the state along the vector field k-^. Now, the 
1-parameter family of diffeomorphisms generated by the vector field d/d/i has a well-defined 
action also on the Hilbert space Ti^in^ under consideration. Therefore, Eq ()2.15j) holds 
also for states ^{fi) in H^j^ provided exp kd/dfi is interpreted simply as the operator that 
drags the state a unit afiine parameter distance along the vector field k-^. This suggests 
that it is natural to set 

d 



e'2\iJ{jj) = e^^^{fi) (2.16) 
so the right side would be the image of \^ under the finite diffeomorphism obtained by 



moving a unit affine parameter distance along the integral curve of the vector field ^ = ^4-. 



The affine parameter v along this vector field is given by: 

3 2 
V = K sgn{fi) , where K = — j= (2-17) 

3V3V3 

Although the geometrical meaning of this action of expi(/ic/2) is simple, since f is a rather 
complicated function of /i, its expression in the //-representation is complicated: 

e*^ ^(/i) = ^'(sgn(/i) I /i 1 3), where /i = sgn(/i) + (2.18) 

K 
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However, it is well-defined because v{^) is an invertible (and C^) function of /i. Furthermore, 
for /i 3> 1, the action reduces to the form 

e^f ^(/i) ^ + (2.19) 

familiar from the 'standard' LQG analysis where ft is replaced by a constant fio (see Eq 
fl2.15|) ). However, as is clear from the exact expression ()2.18|) . this form with a 'simple 
displacement' of the argument is highly inaccurate for small fi. 

The complicated action of exp z(/ic/2) just reflects the fact that the variable /i determined 
by eigenvalues of p is not well-adapted to the vector field ^. Let us therefore change the 
basis to \v). This basis is more directly adapted to the volume operator V (associated 
with the elementary cell V): 

W = (^)^^4iH (2.20) 

(It follows from ()2.17j) that, like fi, v is dimensionless.) These kets also constitute an or- 
thonormal basis in 'H^^^: {vi \v2) = 5i,i,i,2- The two bases are of course closely related 
because V and p are just functions of one another. In terms of representations, we can sim- 
ply set \l/(f) = ^{fJ')- In the f-representation, the action of expi(/2c/2) is extremely simple 

e^f ^(w) = ^(w + 1) (2.21) 

Therefore, in what follows we will use the v representation. Many of the expressions obtained 
using the new constraint will then have a form rather similar to those obtained in 0, 0| ■ 
But it is important to keep in mind that while those expressions were written in the n 
representation, the ones in this paper are written in the ^-representation. This difference is 
a key reason why the present Hamiltonian constraint is free of the drawback of the one used 
in the literature so far. 

Note that the operator exp i{flc/2) is unitary. We can extend the definition in the obvious 
manner: By setting 

e''''^^{v) = ^{v + k) (2.22) 

for any constant k, we obtain a unitary representation on Tif^^^ of the 1-parameter group of 
diffeomorphisms generated by the vector field ^. On the classical phase space, the function 
/2c/2 generates a 1-parameter family of (finite) canonical transformations. It is precisely the 
lift to the phase space of the 1-parameter group of diffeomorphisms on the p or fi space, 

generated by the vector field /id/d/x. The action of exp i{fic/2) simply promotes this relation 
to the quantum theory. 

Remark: While the geometric interpretation of the action of expi{fic/2) makes its defi- 
nition quite natural, one might nonetheless ask for the relation between our factor ordering 
and that normally used. We will now argue that the choice we have specified in fact origi- 
nates in ordinary quantum mechanics. 

Consider 1-d Schrodinger quantum mechanics, where states ^{x) are generally taken to 
be functions in L^(M, dx). Consider a function f{x)p on the classical phase space. The 
canonical transformation it generates is again the lift to the phase space of the diffeomor- 
phism generated by the vector field ^ = /(x)-^ on the configuration space (i.e., the x-axis). 
One would expect that this geometric action would be carried over also to quantum theory. 
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Is this in fact the case? In the textbook treatments one does not^ take into account such 
geometric considerations but constructs a self-adjoint operator (fp) := {h/2i) (/ ^ + ^ /) 

using a symmetric ordering. Its action {fp)'^ = (h/i) (/ ^^ + | ^ / ^) does not correspond 
to the Lie-derivative of ^ with respect to the vector field ^, which is captured just in the 
first term. (The second term, | ^ / ^, multiplies ^ by the divergence of the vector field 
^, i.e., by the Lie-derivative of the Lebesgue measure with respect to ^.) Hence, the action 

of the 1-parameter group of unitary transformations U{X) '■= expi{fp) on wave functions 
^(x) is different from that of the diffeomorphisms generated by ^. However, this is simply 
because the Lebesgue measure dx fails to be invariant under the diffeomorphism. 

Let us represent quantum states by densities of weight 1/2 — or, more precisely, half- 
forms — so that the scalar product is given simply by (\l'i, \1'2) = without the need 

of a measure. (Thus, '^i^2dx ~ 'I'i\E'2- For details, see, e.g., On these half forms, 

the action of (fp) is precisely that of a Lie derivative; (fp) \I' = (h/i) £^\E'. Therefore, the 
action of the unitary transformation U (A) is given by: 

U{X)^{v) = ^{v + Xh) 

where v is the affine parameter of the vector field ^ on R. Thus, the initial geometric expec- 
tation is indeed realized if one proceeds with the standard factor ordering from Schrodinger 
quantum mechanics, but represents states by half forms. In the polymer representation the 
measure is invariant under the action of any diffeomorphism on the /i-axis. Hence the addi- 
tional 'divergence term' is unnecessary. In this sense, states in H^^'^ are analogs of densities 
of weight 1/2 — or, half forms — in the Schrodinger representation. 



C. Expression of the constraint operator 



We can now collect the results of the last two sub-sections to obtain the quantum con- 
straint operator starting from the classical expression ()2.13|) . We have: 



a 



grav 



sin(/ic) 



242 sgn(/i) 



sin(/ic) A sin(/ic) 
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/^C\ y . ffIC 

— \/ sm — 



sin (/ic) 
(2.23) 



where, for clarity of visualization, we have suppressed hats over the operators 
sin(/2c/2), cos(/2c/2) and sgn(/i)//2^. Let us focus on the operator A first. Some care is 
needed in its evaluation because, while all other operators in this expression are densely 
defined, sgn(/i) is unambiguously defined only on those states $(f ) whose support excludes 
the point v = 0. However, a straightforward calculation shows that 
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(2.24) 



(The expression in round brackets (containing the trignometric functions and V) is already 
diagonal in the \v) basis. Hence there is no factor ordering problem between this part of the 
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operator and the pre-factor involving Since the right hand side vanishes ai v = 0, 

it is in the domain of sgn(/i), whence A is well defined and given by 

= -^/f + nv). (2.25) 

Thus, \v) is an eigenket of A. Since the eigenvalues of A are real and negative, it is a negative 
definite self-adjoint operator on Ti^in^- The form of Cgrsw in ()2.23|1 immediately implies that 
it is also self-adjoint and negative definite on Ti.^^'^ . Its action is given by: 

C'grav ^{v) = U{v) ^{v + 4) + ^v) ^ (v) + f.{v) ^ {v - 4) (2.26) 

with 

hiv) = i^Vt^I^ + 'I ||. + i|-|. + 3|| 

f-{v) = U{v-A) 

foiv) = -f4v)-f4v) . (2.27) 

Thus, the new gravitational constraint is again a difference operator. However, whereas the 
operator used so far in the literature involves steps which are constant (4/io) in magnitude 
in the eigenvalues of p, the new constraint involves steps which are constant in eigenvalues 
of the volume operator V. In the \fi) basis these steps vary, becoming smaller for large /j,. 
Note also that although fl diverges at f = 0, individual operators entering the constraint 
— and hence the full constraint operator itself — are well-defined on the state |f = 0). 

Finally, to write the complete constraint operator we also need the matter part of the 
constraint. For the massless scalar field, in the classical theory it is given by: 

Cm.tt = SttG \p\~hl (2.28) 

Thus, as usual, the non-trivial part in the passage to quantum theory is the function 
However, as with the co-triad operator ()2.10p . this can be used by the method introduced 
by Thiemann in the full theory 3]. In this quantization, there are two ambiguities 
[3, , labelled by an half integer j and a real number i in the range < £ < 1 . As in P| , 
following the general considerations in 0,113], we will set j = 1/2. For £, a general selection 
criterion is not available and values i = 1/2 and i = 3/4 have been used varyingly in the 
literature. Qualitative features of our results do not depend on this choice. Since i = 1/2 
makes expressions simpler and since we used £ = 3/4 in Q, for simplicity and variety we 
will set £ = 1/2 in this paper. Then, the point of departure is the identity 



\p\ 2 = sgn(p) 



(2.29) 



on the classical phase space. In the LQC literature, this identity is generally used with a 
constant fj,o in place of fi. However, since volume depends only on p, in the Poisson bracket 
only the derivative with respect to c of the holonomy appears. Therefore, the identity 
continues to hold although /i is a function of /i. ^ 



^ Indeed, one could replace it by any other suitably regular function of ^. We will refrain from doing so 
because the function /2 was already determined for us in the gravitational part of the Hamiltonian and 
choosing another function in an ad-hoc manner for the matter part would make the operator less natural. 
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Again in the passage to quantum theory, some care is needed because of the presence of 
sgn(p) on the right side. However, again the image of the operator defined by the expression 
in the square brackets is in the domain of sgn(p), whence the total operator on the right 
side is densely defined: 



6 



1/2 



i^l/3|^|l/3| |^ + l|l/3 



(2.30) 



For |f| 1, the eigenvalue is given by (6/87r7£pj)-'^/^ sgn(/i) l/il^-*^/^ (1 + 0{l/\fi\'^)) ~ 
sgn(p) IpI"^''^, whence the classical behavior is recovered. However, the operator is well- 
behaved on the ket \v = 0); in fact, as is usual in LQC it is an eigenvector and the eigenvalue 
vanishes. 

Since is a well-defined, self-adjoint operator, we can take its cube to obtain the 

action of It is diagonal in the v representation, with action: 

^— — - / c \ 3/2 

N~^^M=hr^ B{v)^{v) (2.31) 



PI 



where 



B(v) 



K \v\ 



+l|l/3_L,_l|l/3 



Collecting these results we can express the total constraint 



C ^{V) = (4,av + C'matt) "^{V) = , 



(2.32) 



(2.33) 



as follows: 



[B{v)]-'^ {C^{v) m{v + 4, 0) + C°{v) m{v, 0) + C-{v) m{v - 4, 0)) 



(2.34) 



(2.35) 



where the coefficients C {v) and C°{y) are given by: 

C+[v) = — - — \v + 2\ \\v + 1\ - \v + 3\\ 
8 

C-{v) = C+{v-4) 

This is the Hamiltonian constraint we will work with in the remainder of the paper. As 
discussed in P|, the form of this constraint is similar to that of a massless Klein-Gordon 
field in a static space-time, with playing the role of time and the difference operator B of 
the spatial Laplace-type operator. Hence, the scalar field can again be used as 'emergent 
time' in the quantum theory. We will examine the operator B in some detail in sections 
IIVI and 13 Finally, in the above construction we made a factor ordering choice, the most 
significant of which is to write the gravitational part of the constraint as sin ficA sin fic by 
splitting the sin^/ic term jilj (see ()2.13j) and ()2.23|l ). Since A is self-adjoint and negative 
definite (and sin /2c is self-adjoint), this ordering directly endows the same properties on 
the gravitational constraint Cgrav However, there are other possibilities which may well be 
better suited in more complicated models. For the model under consideration, we considered 
one other natural candidate in ^ and found that the main results were not affected by the 
change. Therefore we did not use other factor orderings in the numerical simulations with 
the new constraint. 
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III. THE WHEELER-DEWITT THEORY 

In this section we will discuss the WDW limit of LQC in which effects specific to quantum 
geometry in the difference equation ()2.34|1 are ignored. Although this limiting theory is 
straightforward, we present it in some detail because it provides a simpler and more familiar 
setting for constructing the physical Hilbert space, Dirac observables and semi-classical 
states and because this discussion will enable us to compare and contrast the WDW theory 
with LQC in detail. 

The section is divided into two parts. In the first we obtain the WDW limit of ()2.34j) and 
its general solution. In the second we construct the physical sector of the theory, using the 
scalar field as emergent time and show that the big bang singularity is not resolved in the 
WDW limit. 

A. The WDW constraint and its general solution 

To obtain the WDW limit of ()2.34|) we will follow the same procedure that we used in Q . 
In particular, since the geometrodynamical description is insensitive to the choice of triad 
orientation — i.e., to the sign of v — we will restrict ourselves to wave functions '^{v) which 
are symmetric under v — >■ —v. As in jlj, the WDW analogs of the LQC quantities will be 
written with an underbar. 

Let us begin by setting: 

^^i^v) := {\v\ \\v - l\ - \v + {m{v + 2)-m{v-2)) , (3.1) 
so that the gravitational part (j2.26|) of the Hamiltonian constraint can be written as 

^^rav nv) = [X{v + 2) - x{v - 2)] . (3.2) 

To obtain the WDW limit, let us assume |f | ^ 1 and \I'(f) is smooth. Then, we have: 

C.„*(„)^5«-4,;f -i(K,|^f),.0,„-f;i) (3.3, 

where n > 3. Thus, if we restrict ourselves to wave functions \E' which are slowly varying in 
the sense that the second term on the right hand side is negligible compared to the first, we 
obtain the WDW limit of the gravitational part of the constraint: 

^(IH^). (3.4) 

This approximation is not uniform because the terms which are neglected depend on ^. 
We will show at the end that these assumptions are realized in a self-consistent manner on 
semi-classical states of interest. For now we only note that C'^^'^ is self-adjoint and negative 
definite on the kinematic Hilbert space = L^(]R, dv) of the WDW theory. 

Finally, recall that the full constraint is given by 

C = + SvrG (^^^^ j B{v)pI (3.5) 
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where the function B{v) is defined in ()2.32|) . This function is well-approximated by B{v) 
K\v\~^ for |f I 3> 1. Hence, the WDW limit of the full constraint is given by: 



dl^{v,(l)) = 12ttG vd^{vd^^{i 

=: -e^{v,(f)) (3.6) 

The operator commutes with the 'parity' operator 11 which fiips the triad orientation: 
Il^{v) = :^{-v).'^ Hence it preserves the space of wave functions under consideration, 
namely the eigenspace of H with eigenvalue +1. 

Note that if, as is in geometrodynamics, we were to represent quantum states as functions 
^{a, (p), the operator 9 would have the action: 9;^ = (47rG'/3) {ad / da) {ad/ da) {'^). Now, in 
geometrodynamics, the classical constraint has the form G^^paPb = where G^^ is the De- 
Witt metric on the 2-dimensional configuration space spanned by the pair (a, 0) . In quantum 
theory, the natural choice of factor ordering is to use the Laplace-Beltrami operator of G^^ 
23- Then, the WDW equation has the form: & = (47rG/3) {ad /da) {ad /da)^ lil, which 



is precisely our equation ()3.(i|l . Thus the WDW limit of the current, 'improved' Hamilto- 
nian constraint, automatically yields the 'natural' factor ordering of the WDW theory. This 
provides an indirect support for the factor ordering in LQC that led us to ()2.34p . 

Since the form of the WDW limit is symmetric in and v we could use either the volume 
of the universe or the scalar field as the 'emergent time' and the other variable as the true 
dynamical degree of freedom. However, as noted in (^=1) closed model, is 

better suited as the time variable. More importantly, the form of the LQC constraint ()2.34|) 
is such that it would be technically difficult to regard v as time beyond the WDW limit. 
Therefore, we will regard as time so that the dynamics is described by the evolution of 
the volume of the universe with respect to this emergent time, 0. As emphasized in , it is 
not essential to make a specific choice of an 'emergent time'; one can construct the physical 
sector of the theory without making a choice. However, choosing as time provides a 
heuristic understanding of the intermediate stages of the procedure and lets us interpret the 
final results in the language of the more familiar notion of 'evolution' rather than in terms 
of a 'frozen formalism'. 

The form of the WDW equation ()3.6p is very simple; if we set x = \nv, the 
WDW constraint assumes the form of a massless Klein Gordon equation in the fiat space 
spanned by x and 0. However, to bring out similarities and contrasts between this 
WDW theory and the LQC description presented in sections |TV| and El we will use the 
same steps as the ones used there. 

Let us first note that the operator 9 is positive definite and self-adjoint on the Hilbert 
space L'^{W, K~^B_{v)dv) where the subscript s denotes the restriction to the symmetric 
eigenspace of H. Its eigenf unctions ej^ with eigenvalue ci;^(> 0) are 2-fold degenerate on this 
Hilbert space. Therefore, they can be labelled by a real number k: 

e,{v):=^e^''^\^\ (3.7) 



In the full classical theory, the orientation reversal is a phase space symmetry because both triad e° and 
its conjugate variable fc^ :— kate^^ change signs under this transformation. In the present model, this cor- 
responds to the symplectomorphism {c,p) (— c, —p) which is then naturally lifted to an automorphism 
of the quantum algebra. The mapping 11 is the induced action of this automorphism on states. 
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where k is related to tu via u = \/l2TTG\k\. On L^l-'^) K~^B_(v)dv), they form an orthonormal 
basis: (e^, e'^) = 5{k, k'). A general solution to ()3.6|) with initial data in the Schwartz space 
of rapidly decreasing functions can be written as 

/oo 
dk ^+{k) efc(t;)e*"^ + ^_ (k) ei^{v)e-''^'f' (3.8) 
-oo 

for some ^±{k) also in the Schwartz space. Following the terminology generally used in the 
Klein-Gordon theory, the solution will be said to be 'incoming' ('contracting') if \^±(fc) have 
support only on the positive half of the k axis and 'outgoing' or ('expanding') if they have 
support only on the negative half of the k axis. If \&_(/c) vanishes, the solution will be said 
to be of positive frequency and if '^+{k) vanishes it will be said to be of negative frequency. 

As usual, the positive and negative frequency solutions satisfy first order 'evolution' 
equations, obtained by taking a square-root of the constraint ()3.6|) : 

T^^M^,^) = ^/Q:^{v,4>). (3.9) 

If f{v) is the initial data for these equations at 'time (p = 4'o\ the solutions are given by: 

^^{v, 0) = e^'V^ f{v) (3.10) 



B. Physics of the WDW theory 

Solutions ()3.8j) to the WDW equation are not normalizable in Ti^^ (because zero is in 
the continuous part of the spectrum of the WDW operator). Our first task is to endow 
the space of these physical states with a Hilbert space structure. There are several possible 
avenues. As in ^ , we will beg in with one that is somewhat heuristic but has direct physical 
motivation. The idea [i^, 125^ is to introduce operators corresponding to a complete set of 
Dirac observables and select the required inner product by demanding that they be self- 
adjoint. In the classical theory, since is a constant of motion, it is a Dirac observable. 
While V is not a constant of motion, on each dynamical trajectory f (0) is a monotonic 
function of (f), whence v \(f,=(j,^ is a Dirac observable for any fixed 0o- These form a complete set 
and it is straightforward to write the corresponding quantum operators. However, because 
we are interested only in states '${v,(f)) which are symmetric under v —v, it suffices to 
consider \v\,f,^ in place of f |</,„. Now, since commutes with the WDW operator in ()3.6p . 
given a (symmetric) solution ^!^(f , 0) to ()3.6|) . 




(3.11) 



is again a (symmetric) solution. So, we can just retain this definition of p^ from 
The Schrodinger type evolutions ()3.10|) enable us to define the other Dirac observable \v\(i,^'. 
Given a (symmetric) solution ^(f , 0) to ()3.6p . we can first decompose it into positive and 
negative frequency parts ^±{v,(f)), freeze them at = 0o, multiply this 'initial datum' by 
\v\ and evolve via ()3.10|) : 

\v\4>o'^{v, 0) = e*v^('^-'^°) l^^|:*+(^^, 0o) + e-^v^ \v\:^_{v, 0^) (3.12) 
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The result is again a (symmetric) solution to the WDW equation ()3.6p . Furthermore, both 
these operators preserve the positive and negative frequency subspaces. Since they constitute 
a complete family of Dirac observables, we have superselection. In quantum theory we can 
restrict ourselves to one superselected sector. In what follows, for definiteness we will focus 
on the positive frequency sector and, from now on, drop the suffix +. 

We now seek an inner product on the space of positive frequency solutions ^(f , (p) to 
(|3.1Up (invariant under the v reflection) which makes p^ and l^l,^^ self-adjoint. Each of these 
solutions is completely determined by its initial datum ^{v, (j)o) and the Dirac observables 
have the following action on the datum: 

\vUoM'"^<P) = \'"\M'"Ao). and P4>MvAo) = fi^^{v,<i)o) . (3.13) 

Therefore, it follows that (modulo an overall rescaling,) the unique inner product which will 
make these operators self-adjoint is just: 

(^ill2)phy= / dvB{v)Mv.(t>)Mv.(t>) (3.14) 
J 4>=(l>o 

(see e.g. jl^, HH). Note that the inner product is conserved, i.e., is independent of the 
choice of the 'instant' = 0o- Thus, the physical Hilbert space Tiph™ ^■^ the space of positive 
frequency wave functions ^{v, 0) which are symmetric under the v reflection and have a 
finite norm, defined by ^3.14\ l- The procedure has already provided us with a representation 
of our complete set of Dirac observables on this 'Hp^^: 

I^^Uo = e*v^^'^"'^°^ |w|:*(t',0o), and Pc^^{v,(f)) = h^/e^{v,(j)) . (3.15) 

Arguments of 0| can be repeated to show that the same representation of the algebra of Dirac 
observables can be obtained by the group averaging method ^26^ ,27,] which is mathematically 
more stream-lined. 

Finally, using the physical Hilbert space and this complete set of Dirac observables we 
can now introduce semi-classical states and study their evolution. Let us fix an 'instant of 
time' = 0o and construct a semi-classical state which is peaked at p^p = and f = v*. 
Since we would like the peak to be at a point that represents a large classical universe, we 
are led to choose f * ^ 1 and (in the natural classical units c=G=l) p*^ ^ h. In the closed 
[k = 1) models for example, the second condition is necessary to ensure that the universe 
expands out to a size much larger than the Planck scale. At 'time' = 0o, consider the 
state ^ 

^(t;,0„) = I dk ^<{k) efc(f) e''^^'^°-'^*\ where ^!{k) = e'^^^ . (3.16) 



Here k* = —p*^/V12TiGh? and 0* = —\n\v*\/^/TJT2nG + 0o- It is easy to evaluate the 

integral in the approximation uj = — VlSvrG k — which is justified because \E'(fc) is sharply 
peaked at k* and k* <^ — 1 — and calculate mean values of the Dirac observables and their 
fluctuations. One finds that, as required, at = 0o the state is sharply peaked at values 
v*,p^. The above construction is closely related to that of coherent states in non-relativistic 
quantum mechanics. The main difference is that the observables of interest are not v and 
its conjugate momentum but rather v and p^ — the momentum conjugate to 'time', i.e., the 
analog of the Hamiltonian in non-relativistic quantum mechanics. 
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We can now ask for the evolution of this state. Does it remain peaked at the classical 
trajectory defined by = and passing through v = v* at (j) = (j)o'^ This question is 
easy to answer because ()3.10p implies that the (positive frequency) solution to ^{v, 0) ()3.6|) 
defined by the initial data (j3.16p is obtained simply by replacing 0o by in (j3.16|) ! Since a, 
the measure of dispersion in ()3.1(ij) . does not depend on 0, it follows that ^{v, 0) continues 
to be peaked at a trajectory 

which is precisely the classical solution of interest. This is precisely what one would hope 
during the epoch in which the universe is large. However, the property holds also in the 
Planck regime and, in the backward evolution, the semi-classical state simply follows the 
classical trajectory into the big-bang singularity. (Had we worked with positive k*, we would 
have obtained a contracting solution and then the forward evolution would have followed 
the classical trajectory into the big-crunch singularity.) In this sense, the WDW evolution 
does not resolve the classical singularity. 

We will show in sections ll VI and M that the situation is very different in LQC This can 
occur because the WDW equation is a good approximation to the discrete equation only for 
large v. Furthermore, as noted before, the approximation is not uniform but depends on the 
state: in arriving at the WDW equation from LQC we had to neglect \l/ dependent terms 
of the form 0{v"'~^^^) for n > 3. For semi-classical states considered above, this implies 
that the approximation is excellent for v ^ k* but becomes inadequate when the peak of 
the wave function lies at a value of v comparable to k*. Then, the LQC evolution departs 
sharply from the WDW evolution. We will find that, rather than following the classical 
trajectory into the big bang singularity, the peak now exhibits a bounce. Since large values 
of k* are classically preferred, the value of v at the bounce can be quite large. However, as 
remarked in Sec. HJ we will find that the matter density at the bounce point is comparable 
to the Planck density, independent of the value of k*. 

Remark: In the above discussion for simplicity we restricted ourselves to eigenfunctions 
ek{v) which are symmetric under v ~v from the beginning. Had we dropped this 
requirement, we would have found that there is a 4-fold (rather than 2- fold) degeneracy in 
the eigenfunctions of 6. Indeed, if 9{v) is the step function (6'(f ) = if f < and = 1 if 
V > 0), then 9{v)e\^, 0{v)e_^^, ^{~'v)§.\k\y continuous functions of v which 

satisfy the eigenvalue equation (in the distributional sense) with eigenvalue uj'^ = 12nGk^. 
This fact will be relevant in the next section. 

IV. ANALYTICAL ISSUES IN LOOP QUANTUM COSMOLOGY 

We will now analyze the model using LQC. Since the form of the LQC 'evolution equation' 
is very similar to that of the WDW theory, we will be able to construct the physical Hilbert 
space and Dirac observables following the ideas introduced in section ITll Bl 
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A. Emergent time and the general solution to the LQC Hamiltonian constraint 

Recall that the LQG Hamiltonian constraint is given by Eq ()2.34|) : 



where the coefficients C^,C° are given by ()2.H5|1 .^ Since the operator B acts only on the 
argument v of '^{v^cj)), we have a neat separation of variables. The form of the constraint 
now suggests that it is natural to regard (f) as emergent time. To implement this idea, let 
us introduce an appropriate kinematical Hilbert space for both geometry and the scalar 
field: V}^^^ := L2(RB„hr, 5(v)d/iBohr) ® ^^(R, d0). Since is to be thought of as 'time' and 
V as the genuine, physical degree of freedom which evolves with respect to this 'time', we 
chose the standard Schrodinger representation for but the 'polymer representation' for v 
to correctly incorporate the quantum geometry effects. This is a conservative approach in 
that the results will directly reveal the manifestations of quantum geometry. Had we chosen 
a non-standard representation for the scalar field, these effects would have been mixed 
with those arising from an unusual representation of 'time evolution' and, furthermore, 
comparison with the WDW theory would have become more complicated. (However, the 
use of a 'polymer representation' for may become necessary to treat inhomogeneities in 
an adequate fashion.) 

The form of ()2.34|) is the same as that of the WDW constraint ()3.6|) and G is again a 
positive self-adjoint operator, now on L^(MBohr, -B(^^)d/iBohr)- The main difference is that 
while the WDW G is a differential operator, the LQC G is a difference operator. This gives 
rise to certain technically important distinctions. For, now the space of physical states — i.e. 
of appropriate solutions to the constraint equation — is naturally divided into sectors each of 
which is preserved by the 'evolution' and by the action of our Dirac observables. Thus, there 
is super-selection. Let C\e\ denote the 'lattice' of points {|£:| +4n, n E Z} on the v-axis, jC-\e\ 
the 'lattice' of points { — 1^1 + 4?t,, n E Z} and let Cs = '^lei U -^^-lei where as usual Z denotes 
the set of integers. Let Ti.^^^^ , Ti.^^^^ and Hf^'' denote the subspaces of L^(RBohr, )d/iBohr) 
with states whose support is restricted to lattices C\e\,C-\e\ and Ce- Each of these three 
subspaces is mapped to itself by G which is self-adjoint and positive definite on all three 
Hilbert spaces. 

However, since 'H^^^^ and Ti.^^^^ are mapped to each other by the parity operator H, only 
-^grav jg invariant by H. Now, because H reverses the triad orientation, it represents 
a large gauge transformation. In gauge theories, we have to restrict ourselves to sectors, 
each consisting of an eigenspace of the group of large gauge transformations. (In QCD in 
particular this leads to the 6 sectors.) The group generated by H is just Z2, whence there 
are only two eigenspaces, with eigenvalues ±1. Since there are no fermions in our theory, 
there are no parity violating processes whence we are led to choose the symmetric sector 
with eigenvalue +1. Thus, we are primarily interested in the symmetric subspace of Tif^^; 
the other two Hilbert spaces will be useful only in the intermediate stages of our discussion. 

Our ffist task is to explore properties of the operator G. Since it is self-adjoint and 
positive definite, its spectrum is non-negative. Therefore, as in the WDW theory, we will 



^ Note that this fundamental evolution equation makes no reference to the Barbero-Immirzi parameter 7 
or the area gap A. This is the equation used in numerical simulations. To interpret the results in terms 
of scale factor, however, values of 7 and A become relevant. 



9j^(t;,0) = [B{v)Y^ {C^{v)^{v + A,(t)) + C\v)^{v,(t)) + C-{v)^{v-A,(t))) 




(4.1) 



18 



denote its eigenvalues by uP' . Let us first consider a generic e, i.e., not equal to or 2, so 
that the lattices C.±\e\ are distinct. On each of the two Hilbert spaces T^^^-, we can solve for 
the eigenvalue equation Q e^{y) = uj"^ e^{v), i.e., 

C+iv)e^iv + 4) + C"iv)eM + C-{v)eUv - 4) = uj^B{v)eM (4-2) 

Since this equation has the form of a recursion relation and since the coefficients C^{v) never 
vanish on the 'lattices' under consideration, it follows that we will obtain an eigenfunction 
by freely specifying, say, \E'(t;*) and \E'(t>* + 4) for any v* on the lattice C\s\ or C-.\s\. Hence 
the eigenf unctions are 2-fold degenerate on each of and Ti.^^^^. On TCf^^, therefore, 

the eigenf unctions are 4-fold degenerate as in the WDW theory. Consequently Ti.f^'^ admits 
an orthonormal basis ef, where the degeneracy index / ranges from 1 to 4, such that 

(e^le^') = 6j,j>6{uj,u;'). (4.3) 

(The Hilbert space TCf^^ is separable and the spectrum is equipped with the standard 
topology of the real line. Therefore we have the Dirac distribution 6{uj, uj') rather than the 
Kronecker delta 5^^^^^l.) As usual, every element \E'(f) of Tif^"'' can be expanded as: 

-^{v)= j duj^>i{uj)el{y) where ^ i{uj) = (e^|^) , (4.4) 
Jo 

The numerical analysis of section |3 and comparison with the WDW theory are facilitated 
by making a convenient choice of this basis in Tif^^, i.e., by picking specific vectors from 
each 4 dimensional eigenspace spanned by e^(f). To do so, note first that, as one might 
expect, every eigenvector ef,(t') has the property that it approaches specific eigenvectors of 
the WDW differential operator 6 as f ^ ±oo. (The precise rate of approach is discussed 
in section IVAl ) The idea is to use this fact to synchronize the basis in LQC with that in 
the WDW theory. However, in general the two limiting WDW eigenfunctions are distinct. 
Indeed, because of the nature of the WDW operator B, its eigenvectors can be chosen 
to vanish on the entire negative (or positive) f-axis; their behavior on the two half lines is 
uncorrelated. (See the remark at the end of section UllBl ) Eigenvectors of the LQC operator 
O on the other hand are rigid; their values at any two lattice points determine their values 
on the entire lattice jC±\s\. Therefore, the synchronization can be done only at one end and 
we choose to do so at f = oo. As in the WDW theory, let us introduce a real variable k 
satisfying u"^ = (127rG)/c^ and use k in place of u to label the orthonormal basis. Then the 
basis of interest will be the following: 

i) Denote by e^|^|(f) the basis vector in Ti-^^J^ with eigenvalue cj^, which is pro- 
portional to the WDW e_|^|(f) in the limit f ^ oo; (i.e., it has only 'outgoing' 
or 'expanding' component in this limit); 

ii) Denote by ej^| (v) the basis vector in 'H^^^^ with eigenvalue a;^ which is orthog- 
onal to e^|^,|(f). (Since eigenvectors are 2-fold degenerate in each of T^^\'J\, the 
vector ej^|(f) is uniquely determined up to a multiplicative phase factor.) 

We thus have an orthonormal basis in Hf^^ with A; G M: (e^|e^/) = 5{k, k'), and (e^|e^,) = 
0. The four eigenvectors with eigenvalue uj"^ are now e|^|,e^|j.| which have support on the 

'lattice' and e|^|, e~|^| which have support on the 'lattice' As we will see in section 

IV Al this basis is well-suited for numerical analysis. 
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Since physical states use only the symmetric sector of Tie, in this analytic discussion we 
will be interested only in the symmetric combinations 

= \ (etiv) + eti-v) + e-,{v) + e-,{-v)) (4.5) 

of the basis vectors which are invariant under 11. Any symmetric element \E'(f ) of Tif^'^ can 
be expanded as 

POO 

<i/{v) = dk ^{k) ei"\v) (4.6) 



We can now write down the general symmetric solution to the quantum constraint ()2.34|) 
with initial data in Hf^^ : 

/OD 
dk [^+{k)ei'\v)e''^'^ + ^.{k)ei'\v)e~'^'^] (4.7) 
-oo 

where \E'±(/c) are in L^(M, rf/c). As v — *• ±oo, these approach solutions ()3.8|) to the 
WDW equation. However, the approach is not uniform in the Hilbert space but varies 
from solution to solution. As mentioned in section llll Bl the LQC solutions to ()2.34j) which 
are semi-classical at late times can start departing from the WDW solutions for relatively 
large values of v (more precisely, at f ~ A;). 

As in the WDW theory, if \E'_(fc) vanishes, we will say that the solution is of positive 
frequency and if ^'+(fc) vanishes we will say it is of negative frequency. Thus, every solution 
to ()2.34|) admits a natural positive and negative frequency decomposition. The positive 
(respectively negative) frequency solutions satisfy a Schrodinger type first order differential 
equation in 0: 

= v^^± (4.8) 

but with a Hamiltonian (which is non-local in v). Therefore the solutions with initial 
datum \E'(t>,0o) = f±{v) are given by: 

^±(^^,0) = e±^^(^-*°)/±(t;,0) (4.9) 

So far, we considered a generic e. We will conclude by summarizing the situation in the 
special cases, e = and e = 2. In these cases, differences arise because the individual lattices 
are invariant under the reflection v — > —v, i.e., the lattices C\e\ and coincide. As before, 
there is a 2-fold degeneracy in the eigenvectors of on any one lattice. For concreteness, let 
us label the Hilbert spaces and choose the basis vectors e^(f), with /c G M as above. 

Now, symmetrization can be performed on each of these Hilbert spaces by itself. So, we 
have: 

ei'\v) = ^ietiv)+ eti-v)) (4.10) 

However, the vector e|^|(f) coincides with the vector e^|^^|(f) so there is only one symmet- 
ric eigenvector per eigenvalue. This is not surprising: the original degeneracy was 2-fold 
(rather than 4-fold) and so there is one symmetric and one anti-symmetric eigenvector per 
eigenvalue. Nonetheless, it is worth noting that there is a precise sense in which the Hilbert 
space of symmetric states is only 'half as big' in these exceptional cases as they are for a 
generic e. 
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For e = 2, there is a further subtlety because vanishes at f = —2 and C~ vanishes at 

V = 2. Thus, in this case, as in the WDW theory, there is a decouphng and the knowledge of 
the eigenfunction e^(f ) on the positive f-axis does not suffice to determine it on the negative 

V axis and vice-versa. However, the degeneracy of the eigenvectors does not increase but 
remains 2-fold because the ()2.H4|1 now introduces two new constraints: C^(±2)e^(±6) = 
[uj^B{±2) — C°(±2)]e^(±2) = 0. Conceptually, this difference is not significant; there is 
again a single symmetric eigenfunction for each eigenvalue. 



B. The Physical sector 

Results of section IIV Al show that while the LQC operator B differs from the 
WDW operator B in interesting ways, the structural form of the two Hamiltonian constraint 
equations is the same. Therefore, apart from the issue of superselection sectors, introduc- 
tion of the Dirac observables and determination of the inner product either by demanding 
that the Dirac observables be self-adjoint or by carrying out group averaging is completely 
analogous in the two cases. Therefore, we will not repeat the discussion of section lIllBI but 
only summarize the final structure. 

The sector of the physical Hilbert space Ti-phy labelled by e G [0, 2] consists of positive 
frequency solutions to (|4.8p with initial data \l/(u,0o) in the symmetric sector of 

Hgrav- Eq. ()4.7|) implies that they have the explicit expression in terms of our eigenvectors 
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*(i;,0)=/ dk'^{k)el''{v)e''^^ , (4.11) 



where, as before, ui'^ = 12'nGk'^ and e^''^(f) is given by (j4.5|) and (j4.1(J|) . By choosing appro- 
priate functions '^{k), this expression will be evaluated in section IV Al using Fast Fourier 
Transforms. The resulting \l/(f , 0) will provide, numerically, quantum states which are semi- 
classical for large v. The physical inner product is given by 

(^i|^2). = B{v)^,{v,ct>o)^2{v,ct>o) (4.12) 

t)e{±|£|+4n;nGZ} 

for any The action of the Dirac observables is independent of e, and has the same form 
as in the WDW theory: 

|{)|<^„^(t;,0) = e^-^(^-^°)|t;|^(t;,0o), and p^^(t;, 0) = - ^^Ml . (4.13) 

The kinematical Hilbert space is non-separable but, because of super-selection, each 

physical sector 'Hphy is separable. Eigenvalues of the Dirac observable |-0|^^ constitute a 
discrete subset of the real line in each sector. The set of these eigenvalues in distinct 
sectors is distinct. Therefore which sector actually occurs is a question that can be in 
principle answered experimentally, provided one has access to microscopic measurements 
which can distinguish between values of the scale factor which differ by ~ 1.3£pi. This will 
not be feasible in the foreseeable future. Of greater practical interest are the coarse-grained 
measurements, where the coarse graining occurs at significantly greater scales. For these 
measurements, different sectors would be indistinguishable and one could work with any 
one. 
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Remark: The procedure used in this section is quite general in the sense that it is apph- 
cable for a large class of systems. In this sense, the physical Hilbert spaces Ti.^Yiy constructed 
here are natural. However, using the special structures available in this model, one can also 
construct an inequivalent representation which is closer to that used in the WDW theory. 
The main results on the bounce also hold in that representation. See Appendix C of 0. 



V. LQC: NUMERICAL ISSUES 

In this section, we will find physical semi-classical states in LQC and analyze their proper- 
ties numerically. This section is divided into three parts. In the first we study eigenf unctions 
e^(/i) of and introduce a method of constructing a 'general' physical state by direct eval- 
uation of the right side of 1)4.1111 . In the second part we solve the initial value problem 
starting from initial data at = 0o, thereby obtaining a 'general' solution to the difference 
equation 1)4.11) . In the third we summarize the main results. Readers who are not interested 
in the details of simulations can go directly to the third subsection. 

A large number of simulations were performed within each of the approaches by varying 
the parameters in the initial data and working with different lattices Cg introduced in Sec. 
IIVI They show the robustness of final results. To avoid making the paper excessively long, 
we will only show illustrative plots. 



A. Eigenfunctions of the Q operator, direct evaluation of the integral 

representation of the state 

Here we will establish properties of the general eigenfunctions e^(f) of introduced 
in section IIV Al briefly present the method of explicit calculation of eigenfunctions in the 
symmetric sector and construct semiclassical states through direct evaluation of their integral 
representation given by ()4.7)) . This method is an exact analog of the one used in [l| and plays 
only an auxiliary role, namely to incorporate the lattice £e=o which is difficult to handle by 
the evolution method described in Sec. IV Bl Therefore our discussion will be brief. 

Let us choose a lattice, say C\e\. Consider an eigenfunction e^^(f ) of C^^^ , i.e. a solution 

to Qe^Jv) = co'^e^(w). Since the left side of this equation approaches [C^Ji^^ )e^\c^^^ as 
|f I — > cxD one expects every solution e^ji^v) to ()4.2|) to converge to some WDW eigenfunction 
e^(f) in this limit. This expectation was verified numerically, following the same procedure 
as the one used in |l|. Numerical simulations also provided the rate of approach^: for 
any e^{v) there exist eigenfunctions e^_+(f), e^_(^^) of (corresponding to the same 

eigenvalue cj^) such that: 




+ 0{^) , foTV>0, 

+ 0{^) , foTv<0. ^ ■ ' 



An example of e^{v) and its approach to e^ _^_{v) for large positive v is presented in fig. [T] 



The numerical tests have shown that the quantity v'^\ei^{v) — ±(w)| is bounded. 
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FIG. 1: Eigenfunction euj{v) of (denoted by +) for e = and lij = 20 is compared with eigenfunc- 
tion e^{v) of the operator in WDW theory (sohd hne). At large v both eigenfunctions approach 
each other, however they differ at small v. 

Each eigenfunction e^(f ) can be expressed as the combination of basis functions e^(f ) 
defined via Eq. (|3.7p : 

eUv) ^ Ae^.^iv) + i?e_|,|(t;), eM C e^.^iv) + De_^,^iv) ^ 

Since the eigenfunctions e^(f ) of 6 are determined on the entire lattice C\s\ by their values 
on (at most) two points, the WDW limits for positive and negative v are not independent. 
Thus, the coefficients C, D are uniquely determined by values of A, B (and vice versa). 

The property (j5.1|) allows us to construct a basis e^^]{v) for the physical (symmetric) 
sector following step by step the procedure presented in [l|: 

1. One first constructs the basis functions e^^^{v) (supported on C±\e\ respectively) by 
solving ()4.2|) on domains C±\e\ fl [— for some fixed, large, positive G C±\£\. 
The initial conditions are specified by demanding that the values of e^|^|(f) agree 
with those of e^(f) at f = + 4. As in the case studied in eigenfunctions 
are enormously amplified on the negative v side and their limits for w -C — 1 are the 
combinations Ce]^^ + -De_|^| of WDW basis functions with coefficients of almost equal 
absolute value i.e., \C^\ ~ 

2. Next, the functions e^^^{v) are used to calculate symmetric eigenfunctions e^|^^|(i;) 
via (j4.5j) . Again their behavior is dominated by properties of e^|^|(w) for f < 0. In 
particular their limit for large v is composed of incoming {k > 0) and outgoing {k < 0) 
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WDW eigenfunctions. The coefficients of decomposition with respect to the basis 
defined by ()3.7p have almost equal absolute value for e ^ 0,2: 

e_|fc| ^z^{e e\k\+e e_|fc|) , (5.3) 

where are some complex constants satisfying l^+l \z~\, while the phases are 
functions of e and uj. For e = and e = 2, = \z~\ and = a". 

These numerically constructed symmetric eigenfunctions can be now used to build the 
semiclassical states, using ()4.1H) . Our aim is to construct a state sharply peaked at a phase 
space point lying on a classical trajectory of the expanding universe at late times. The form 
of the physical inner product and the expression of general physical state given by ()4.1ip 
provides a natural choice of 

^!{k) := e-^^^e-*^-^* (5.4) 

where a is a suitable small spread. The resulting physical state \l/(f,0) is peaked at = 

= —yVlTiGWk. The value of the Dirac observable |f 1,^^ at which the state will be 
peaked at = 0o is determined by 0*. To obtain a semiclassical state at late times we have 
to choose large value of : p'^ h (in units c=G=l). Thus we need to select k* ^ — 1. 

The resulting \1/(A;) has such a small amplitude for A; > that without loss of numerical 
precision it can there be set to zero. Therefore the explicit form of e|*^| is not needed. 
Using expression ()5.4|) in ()4.11|) we obtain final form of the wave function: 

^{v,(f)) = / dke-^^ ei'\v)e''^^'''>'^^-^*K (5.5) 
Jo 

In order to calculate the integral we select a set of k's uniformly distributed within the 
interval [k* — 10a, k* + lOcr], compute the eigenfunctions e'^_^i^^{v) for all k in this set and 
in the chosen finite domain in v, and finally evaluate the integral ()5.5|) using Fast Fourier 
Transform. Numerical simulations were performed with various values of ranging between 
100 and 1000. Typical number of points constituting the interval k ranged from 2048 to 
4096. 



B. Evolution 

As we discussed in Sec. IIV[ the quantum constraint equation can be cast in the form 
of a Klein-Gordon equation in a static spacetime with playing the role of time. It is 
then natural to consider the evolution in terms of the initial value problem in 0. Since 
is a discrete operator in v, the quantum evolution reduces to solving a system of coupled, 
ordinary, 2nd order differential equations. 

To solve the quantum constraint and determine the evolution of the state \l/(f , 0) numer- 
ically, we have to first specify the domain of integration. It was chosen to be |f — e| < 4iV 
where 1 <^ G N is an integer, such that the initial (and, as the evolution shows, also 
subsequent) wave function is negligibly small at the boundary. Nonetheless, since this do- 
main is finite we have specify boundary conditions. We recall that for 1^1 3> 1 the difference 
equation can be very well approximated by the WDW equation. In particular the funda- 
mental equation id^"^ = \/Q^ can be approximated by d^p'^ = (\/l27rG)f(?^\l/. To ensure 
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deterministic evolution, the wave packets , 0) were required to leave (rather than enter) 
the domain of integration. The resulting equation at the boundary was again discretized, 
taking the form: 

= s^?,^G/A{\v\ -2) {^{v.cf)) - ^(t;-4sgn(t;),</))) . (5.6) 

Returning to the problem of evolution on the entire domain the initial data '^\4i^i d^'^\(fi^ 
sX (j) = (po were chosen to be the same as those of the WDW Gaussian semiclassical state 
()3.16|) corresponding to a large expanding universe, peaked at = and f = ^ 1. 
Note that eigenf unctions of the quantum constraint e'^_^^{v) approach a linear combination 
of eigenfunctions of WDW constraint for \v\ 3> 1 fEq. ()5.3| )). Thus in order to keep the 
properties of the initial state specified at |f | ^ 1 as close as possible to those of the solution 
\E'(f , 0) of the 9 operator at large values of |f |, each of the basis functions e_|^| was rotated 
by a phase e~*"±*^'^'^). The phase was found numerically by method analogous to the one 
used in jlf. After its values were determined for variety of different fc, the function of the 
form 

a^{k) = A \n{Bk + C) + D (5.7) 

(where A, B, C, D are real constants) was fitted to them. As the 0th and 1st order of its 
expansion in k corresponds respectively to a constant phase and the shift of origin in these 
terms were subtracted from a"*". The resulting WDW solution then took the form: 

^(t;,0) = j dA;e-^^#e-*°^'e_|;,|(T;)e'"('=)(*-'^*) , (5.8) 

where 

a+'(A;) := a+(A;) -a+(0) - A;9fea+(A;)|o , <f := (po - ^/l/l27iG\n\v*\ . (jSH) 

The right side of ()5.8|) and its derivative w.r.t were calculated numerically to obtain the 
desired initial data. 

This initial data was then evolved using fourth order adaptive Runge-Kutta method. 
To estimate the numerical errors due to discretization in the profiles calculated for 
different step sizes were compared. As the measure of the 'distance' between them we used 
the norm 

11/11(0) := sup |/(t;,,0)| . (5.10) 

\v,-e\<N 

The step sizes were refined until the distance between the results \E'a</,, '^A<f>/2 of integration, 
with step size A0 and A0/2 respectively, satisfied the inequality 

||^A0 - ^A0/2|| < ||^A</,/2|| (A0) S (5.11) 

for a pre-specified, small 6. An example of behavior of the distance for particular calculation 
is presented on fig. |21 They show that the numerical errors manifest themselves mainly 
in phases. The differences between the absolute values of the wave function profiles are 
approximately one order of magnitude smaller. Thus the expectation values and dispersions 
of observables \v\(j, are determined with much better precision than \Ef itself. 

Numerical simulations were performed using 10 different values of which ranged be- 
tween 10'^ and 2 x 10^. The value of v at which the state peaks was always chosen to be 
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FIG. 2: Error functions ||^'(m) ~ ^ll/ll^ll (upper curve) and — |^|||/||^'|| (lower curve) 

are plotted as a function of time steps. Here ^(^m) refers to final profile of wave function for a 
simulation with M time steps. ^ is the limit of final profile as 1/M — > calculated via polynomial 
extrapolation. In both cases, the evolution began at </> = and the final profile is evaluated at 
6 = -1.35. 



greater than 2.5^^. We chose 8 different e sectors, with e values uniformly distributed in 
the interval between an 2, excluding 0. The values of dispersions a/k* varied from 1.7% 
to 11% depending on the value of p^. 

The resulting wave functions \l/(f , 0) of the exact theory were finally used to calculate the 

expectation values (p^), <y|f|(/)^ of observables defined by ()4.13|) . Using the inner product 

(\E'|\E')g given by (j4.12p . their explicit expressions are given by 



vl/) = (vl/|vi/);i J2 Biv)\v\\^{v,^)\' 



(5.12a) 
(5.12b) 



vec. 



where C^^n '■= {v = ±€ + An; —N < n < N} 
The dispersions 



A\v\ 



pI 



\v\4, 



(5.13a) 
(5.13b) 



were also calculated. 
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FIG. 3: The absolute value of the wave function obtained by evolving an initial data using the 
method of (j) evolution. For clarity of visualization, only the values of |^'| greater than 10~^ are 
shown. Being a physical state, ^ is symmetric under v —v. In this simulation, the parameters 
were: e = 2, = 5 x 10^, and Ap^/p'^ = 0.025. 



C. Results 

Results of various numerical simulations can be summarized as follows: 

(1) The states remain sharply peaked throughout the evolution. Their norms are preserved 
under the '0-evolution', providing a consistency check on the numerics. 

(2) The classical trajectory provides a good approximation to the expectation values of 
Dirac observables when the energy density p of the scalar field is small compared to 
a critical energy density pcrit- However, when this value is reached, the expectation 
values cease to be peaked at the classical solution. Instead of following the classical 
trajectory into the big-bang singularity as in the case of the WDW dynamics, the 
LQC state undergoes a quantum bounce. The numerical value of pcrit was the same in 
all simulations, given by pcrit ~ 0.82ppi = 0.82 /G'^h. A physical understanding of this 
phenomenon can be obtained from the modified Friedmann equations which incorporate 
the leading corrections due to quantum geometry effects (see Appendix IB These 
effective equations also provide an analytical expression pcrit = 3/(167r^7'^G^/i) of the 
critical density whose numerical value agrees with that found in the simulations of the 
exact LQC dynamics. 

(3) On further backward evolution p decreases. When it becomes small compared to pcrit, 
the state again becomes sharply peaked on a classical solution which has the same value 
of (p^), but which is contracting in the future. Thus, the quantum geometry effects 
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FIG. 4: The expectation values (and dispersions) of \v\(j) are plotted for the wave function in Fig. 
121 and compared with expanding and contracting classical trajectories. 



lead to a resolution of the big bang singularity. Furthermore, the state becomes semi- 
classical also in the distant past, and the pre and post big-bang branches are joined by 
a 'quantum bridge' by the deterministic evolution of LQC. So long as the initial state 
at late times is chosen to be semi-classical, this scenario, including the value of pcrit, 
is robust. An example of the result of numerical simulation is shown in Fig. |21 The 
comparison between the classical trajectories and quantum evolution is presented in 
Fig. H 

(4) In order to better understand the behavior of energy density during quantum evolution 
we independently analyzed the evolution of expectation values of the energy density 

operator, defined as = p^/2p^|0. Results were consistent with those reported above. 
We found that in all quantum solutions {p^) is bounded above by pcrit ~ 0.82ppi. An 
example of the behavior of the expectation values of p^ is presented in Fig. El 

(5) In the regime where quantum gravity modifications to the classical dynamics are neg- 
ligible, the relative dispersions of Dirac observables were found to remain constant and 
equal on both sides of the bounce point (see Fig. EJ. Thus both semiclassical branches, 
before and after the bounce, are equally coherent. Close to the bounce point the value 
of relative dispersion of \v\^ decreases reaching minimal value exactly at the bounce 
point. However, this decrease does not imply increased coherence in the large density 
region. Indeed using the property that states remain sharply peaked, we can find A0 
from Av/v and estimate the uncertainty product A(j)Ap^. While the product remains 
roughly constant before and after the bounce, it has a small increase near the bounce 
point. 
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FIG. 5: Plot of the behavior of relative dispersion A|r;|/|t;| as function of (p. The upper and lower 
curves correspond to Ap^p/p^ = 0.025 and Ap^p/p^ = 0.018 respectively. As can be seen A|f is 
asymptotically a constant for both expanding and contracting branches. 



Behavior of the dispersion Ap^ is shown in Fig. IHl Its value grows in the classical 
regime in the expanding phase but becomes very small near the bounce point. After 
the bounce, dispersion again increases, before decreasing in the classical contracting 
phase. This peculiar variation in Ap^ near the quantum bounce can be qualitatively 
understood using effective theory discussed in Appendix B. It has its origin in the 
phenomena of super-inflation — i.e. the phase in which H > where H = a/ a is 
the Hubble rate — which occurs in the regime pcrit > P > Pcrit/2 [i^. However, since 
some of the assumptions that underlie the current derivation of the effective equations 
are generally violated in the Planck regime, this 'unreasonable efficacy' of the effective 
equations remains somewhat mysterious. 

VI. DISCUSSION 

This is the second in a series of detailed papers whose goal is to develop a comprehen- 
sive LQG framework to systematically investigate the physics of the Planck regime near the 
big bang. In |lj] we analyzed the homogeneous isotropic model with a massless scalar field 
by introducing several new techniques to construct the physical sector of the theory. We 
established that, thanks to the quantum geometry effects which distinguish LQC from the 
WDW theory, in the backward evolution of states which are semi-classical at late times, the 
big bang is replaced by a quantum bounce. That investigation was based on the quantum 
Hamiltonian constraint that has been used in the LQC literature for the last three years. 
We were able to analyze the physical sector of the theory in detail and firmly establish that. 
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FIG. 6: Plot of the behavior of {p^f,) (points with error bars) and its absolute dispersion (which for 
clarity is multiplied by 50 and is shown by solid curve). Value oip*. = 10, 000 and l^p^/p*. = 0.018. 



although the quantum evolution dictated by that constraint has several desirable features, 
it also has a serious flaw: quantum effects can dominate and signiflcantlv modify classical 
predictions even when the matter density (and curvatures) are low 0, y]. In this paper 
we showed that a physically motivated modiflcation of the quantum Hamiltonian constraint 
overcomes this weakness while preserving the attractive features of the older evolution. In- 
deed, a key strength of this work is that several features of the present simulations are 
qualitatively similar to those of 0. Finally, in both works, detailed numerical analysis 
of dynamics was restricted to those quantum states which are semi-classical at late times. 
Within the model, these states best represent our physical universe and our discussion of the 
quantum bounce refers only to these states, although evolution is well deflned and unitary 
on the full physical Hilbert space. 

We deliberately followed the same organization as ^ to bring out the (similarities and) 
differences between the two quantum evolutions. Physically the differences are critically 
important but mathematically they are rather subtle. This is just as one would expect 
of an improvement that cures a signiflcant limitation while retaining the strengths of the 
older method. In Appendix A we outline the inclusion of the cosmological constant. Again, 
while the previous quantum evolution shows certain departures from the classical theory 
even when the space-time curvature is low \Q] , the quantum evolution generated by the new 
Hamiltonian constraint is free of this drawback. Now the departures occur only in the deep 
Planck regime near the big bang (or the big crunch) and lead to a replacement of the classical 
singularity by a quantum bounce. 

Models discussed so far are too simple to be physically realistic. However, the methods 
introduced are rather general. In particular, as indicated at the end of section III At the '/2- 
strategy' incorporates certain key features of full LQG. Hence, results obtained in this paper 
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provide concrete indications of how quantum geometry can lead to subtle yet important 
departures from the classical theory by making gravity repulsive in the deep Planck regime. 
Glimpses of new physics that may emerge were presented in the last section of We 
direct the reader to that paper also for a detailed discussion of numerical analysis, physical 
ramifications of results and comparisons with other approaches that also feature emergent 
time and/or quantum bounce. 

Finally, as emphasized in section HI a key limitation of the present work is that the theory 
is not obtained through a systematic symmetry reduction of full LQG. This is inevitable 
at present because our understanding of the full quantum dynamics of LQG is still very 
incomplete. The overall framework is constructed by following strategies developed in full 
LQG. However, an important difference arose in the introduction of the Hamiltonian con- 
straint operator because our mini-superspace reduction is carried out by gauge fixing and 
is therefore not diffeomorphism invariant. As a result, in the last step we had to use some 
physical considerations and borrow the expression of the area gap A from the full theory. 
The strategy in LQG is to begin with simple models and work one's way up by incorporating 
the lessons learned along the way. Improved dynamics discussed in this paper is an inter- 
esting example of such lessons. Gonsiderable work is now in progress, also by several others, 
to extend the present analysis to more complicated models that will include anisotropics as 
well as inhomogeneities. 
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APPENDIX A: NON-ZERO COSMOLOGICAL CONSTANT 

In this appendix we will investigate the dynamics of the universe with a massless scalar 
field and a cosmological constant A. Our goal is only to test whether the quantum bounce 
scenario is robust and if the /i-evolution is free of the undesirable features of the /lo-evolution, 
such as departures from the classical theory at low values of total energy densities 0. 
Therefore the numerical simulations are not as refined as in the main body of the paper. A 
more systematic and detailed analysis will appear elsewhere. 

A negative A leads to a classical recollapse of the universe when total energy density, 
including the 'vacuum energy density' due to A, vanishes. If a quantum bounce exists in 
this model then it can lead to a cyclic model of the universe. A positive A is favored by 
current observations. For completeness we will discuss both the cases and ask whether 
general features of the A = analysis persist. 

The classical Hamiltonian constraint is of the form 

Ca = -^c'^ + 2A\p\l +8itG\p\'IpI . (Al) 

From Hamilton's equations it is easy to see that the momentum p,^ is again a constant of 
motion and the scalar field is a monotonic function of time. Hence it is well suited to play 
the role of emergent time in the quantum theory. 



31 



The quantization procedure and the construction of the physical Hilbert space (and ob- 
servables) for the case A 7^ is completely analogous to that in the model with A = 0. The 
quantum Hamiltonian constraint it gives is similar in its form to ()4.H) 

dl^{v,ct>) = -eA^iv,<P) := - 

where Ga is a self-adjoint operator and K is a constant defined in (|2.17p . For A > it fails 
to be positive definite on T^kin- However, since the operator 9| on the left side of ()A2|) is 
negative definite, only the projection of 6a on its positive eigenspace is relevant for solutions 
to the constraint. We can therefore repeat the procedure of section IIVBI decompose the 
solutions into positive and negative frequency parts, and construct the physical inner product 
and observables. The expectation values of \v\^, and their dispersions at an instant are 
again given by ()5.12|) and ()5.13p respectively, where as the norm of the wave function can 
be calculated via ()4.12|) . 

The next step is to construct states which are sharply peaked at late times and compare 
the behavior of the expectation values of observables with classical trajectories. The partic- 
ular properties of the model differ for A < and A > 0. Therefore they have to be handled 
with different methods and we do so separately in sections lA II (A < 0) and lA 21 (A > 0). 
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1. Negative Cosmological Constant 



The classical equations of motion imply that f (0) satisfies the following differential equa- 
tion: 

{3K) 



127iG + )^v'^ . (A3) 



When the cosmological constant is negative the solution to this equation takes the form 

vi<P) = ^^^ (4^yf)3/2^ l cosh(v/T2^(0 - 0.))-i| . (A4) 



Thus, the universe originates at the big bang singularity (for = —00), expands until the 
recoUapse point (at = 0o) where the total energy density (due to the scalar field and the 
cosmological constant) drops to zero and then contracts, reaching the big crunch singularity 
at = +00. 

To determine the quantum evolution we apply the method specified in section IVBl that 

is: 

(i) We first choose on the v axis the domain fl [—Vb,Vb], where Cs is a lattice defined 
in section HV Al for some e 7^ 0, and Vb ^ 1. 

(ii) Next, we specify the initial data, \l/(f , 0o), d^'^i^v, peaked at a phase space point 
representing an expanding, large classical universe. 

(iii) We evolve this data backward, solving ()A2|) numerically as in section IVBl 

(iv) Finally we calculate the expectation values and their dispersions using ()5.12|) and 
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FIG. 7: The expectation values (and dispersions) of |f 1^ for the model with negative cosmological 
constant are compared with classical trajectories. In this simulation, the parameters were: e = 2, 
= 5 X 10^, and A = -3.6 x lO'^, Ap^/p* = 0.034. 



2 / 2 
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We specified the initial state by choosing a Gaussian in (/i, c) peaked at /i* = K~3{v*) 
(see (I2.17|l ) and c*{fi*,p'^), and its time derivative by using the classical trajectory. (Thus, 
we follow 'method I' of section V.B.2 of 0). As discussed there, while this method is not 
as optimal as the one used in main body of this paper, it has the advantage that it does 
not require the knowledge about properties of the eigenfuncions of Ga. The exact forms of 
these initial profiles were the following: 



^{v,(f))\^^ = |t;|ie- '''"2;'*'' e-^(^("^-^*)^^ (A5a) 



2 



5 'roj 



(A5b) 

where f3 := ATT'ji^JS. 

Because the domain of integration in numerical simulations is compact in v, we have to 
provide the boundary conditions. Unlike for A = however the size of classical universe 
is bounded by a maximum value fmax- Therefore, it suffices to choose Vh such that the 
classical recoUapse point is well within the domain and set the 'mirror boundary conditions': 
\1/ = = 0. In numerical simulations Vb was chosen to satisfy the inequality Vb > 1.2t>max 
to make sure that the recoUapse occurs due to dynamics and is not an artifact of the boundary 
conditions. 

An example of the result of simulation is presented in fig. [7| The state remains sharply 
peaked throughout the evolution. The expectation values follow classical trajectory corre- 
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spending to the expansion phase till the total energy density 

Ptot = + PK, where pA := (87rG')~^A , (A6) 

becomes comparable to Pcrit- At Ptot = Pcrit the universe bounces due to repulsive effects 
of quantum geometry, and follows contracting phase of classical trajectory. The agreement 
with classical trajectory remains good through the classical recollapse and expansion phase 
preceding it, until the energy density becomes large again. Then again a bounce is observed. 
In effect the evolution is periodic: each 'cycle' of classical evolution is connected through 
quantum bridge with previous/next cycle. Both the Big Bang and Big Crunch singularities 
are resolved and replaced by quantum bounces. The behavior of fluctuations of Dirac ob- 
servables is also periodic in the sense that the spread at a given point of the trajectory and 
after full cycle are the same. Thus, in the fully quantum evolution of LQC, semi-classical 
states do not loose coherence in evolution from one cycle to another. 



2. Positive Cosmological Constant 

In the case when cosmological constant is positive, classical trajectories, i.e., solutions to 
(jA3|) . take the form 



= v/l2^ (4^^^2^)3/2 sinh(v/l2;^(0 - (/),))-^| . (AT) 



As for A = we have then two types of trajectories. In one the universe starts from a big 
bang singularity, expands and reaches infinite value of v for finite = 0o- In the other the 
universe contracts from the infinite volume state (at a finite = 0o) and reaches the big 
crunch singularity. 

The equation ()A3|) implies in particular that in the region where energy density of scalar 
field is small with respect to vacuum energy density pa, the speed d^v of a wave packet 
following classical trajectory is proportional to f ^. This feature, and the fact that to specify 
proper boundary conditions for evolution in we need to know an explicit form of the square 
root of the Wheeler-DeWitt limit of the positive part of 0a, makes the application of the 
evolution method numerically difficult. Therefore for the construction of semi classical states 
we used the method of direct evaluation of the integral representation of \l/(f,0) specified 
in section IV^ Thus, 

(i) We first calculate the symmetric eigenfunctions e^(f ) of the |9a| operator. 

(ii) Next, we construct the Gaussian state peaked at some fc* and with spread a (see Eq. 

(iii) Finally, we evaluate the integral (|5.5p using fast Fourier transform and calculate the 
expectation values. 

The exact application of this method requires the normalization of e|(f) with respect 
to the inner product given by (j4.12j) . However, our goal here is only to test robust- 
ness of the quantum bounce in presence of a cosmological constant. Therefore, we will 
not construct semi-classical states which minimize uncertainties but work just with states 
which are reasonably sharply peaked. Then, it is enough to use for basis eigenfunctions 
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•= C'(/c)e|(f) where C{k) changes slowly with k, remaining approximately constant 
within interval [k* — 3a, k* + 3a]. To construct such functions we first observe that, for large 
V, e^(f ) show oscillatory behavior and the absolute values of their extremas decrease as \v\ 
increases and set the 'normalization' by requiring equality of maximas closest to the large 
V* chosen initially. 

An example of the result of simulation is presented in figure |H1 The state indeed is sharply 
peaked. As in the case A = 0, the expectation values of \v\^ again reproduce the picture of 
two semiclassical regions, one contracting and the other expanding, connected by a quantum 
bounce. The total energy density (defined via ()A6|) ) at the bounce point equals pcrit- Unhke 
in the A = case however the universe shrinks from v = oo, bounces and again re-expands 
to f = oo, all within a compact interval of 0. 




2*10^ 4*10^ 6*10^ 8*10^ 1.0*10'^ 1.2*10'* 1.4*10'* 



FIG. 8: The expectation values (and dispersions) of for the model with positive cosmological 
constant are compared with classical trajectories. In this simulation, the parameters were: e = 2, 
= 1.6 X 10^, and A = 4.1 x 10~^. Horizontal dashed lines denote asymptotic values of (j) for 
which — > oo. 



To summarize, in both cases (A > and A < 0): 

1. Classical singularity is replaced by a quantum bounce, which occurs when ptot = 
P(t> + Pa = pent- The critical value of energy density was found numerically to be equal 
to 0.82ppi. Thus, the value of pent is independent of the value of cosmological constant. 
Analysis of the modified Friedmann equation on the lines of Appendix B yields the 
same expression, pcrit = V^/ilGrc'^'j^G'^h), of the critical density. 

2. In the dynamics dictated by the Hamiltonian constraint used so far in literature, a 
deviation from the classical behavior resulting in a recoUapse occurs at low energy 
density scales, when cosmological constant dominates The improved dynamics 
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presented in this work is free of this physically undesirable feature. More generally, 
the new dynamics reproduces the standard Friedmann dynamics when p <^ Pcrit- 



APPENDIX B: SOME CONCEPTUAL ISSUES 
1. Effective dynamics 

An effective description of quantum dynamics can be obtained by applying geometric 
methods to quantum mechanics, where the Hilbert space is treated as an infinite dimensional 
phase space which has a structure of a fiber bundle (see, e.g., jl^ E^). The base space of 
this bundle is the classical phase space and fibers are the states with same expectation values 
of the operators for the corresponding canonically conjugate phase space variables. Using 
coherent states we can attempt to find horizontal sections which are preserved under the 
quantum evolution up to some desired accuracy, and can thus obtain an effective Hamiltonian 
which incorporates the leading quantum corrections to the classical dynamics. 

Let us consider the A = case as in the main body of the paper. Then, for LQC, this 
procedure leads to an effective Hamiltonian with leading order terms as j^: 

Kfr = = -^T^bl^ sin^O) + \b{p)pI . (Bl) 

Here -B(p) denotes the eigenvalue of operator given by (j2.31|) . Modifications to the 

dynamics due to behavior of -B(p) become significant for |f | ~ 1. For |f | ^ 1, -B(p) quickly 
approaches the classical value proportional to l/|f | and the rate of approach is given 

by 

g \ 3/2 ^ / 5 ^ / ^ 



Bip) = i + _ +0 — . (B2) 

Neglecting the higher order quantum corrections to B{p) we can obtain the Hamilton's 
equation for v and 0: 

V = {v, n,s} = — - — — = — \6~ ) ^'''^^''^ "'"''^^''^ ■ ^ ^ 

and 



0={0,Kff}= K^y (B4) 

The modified Friedmann equation can be obtained from the vanishing of the Hamiltonian 
constraint ()B1|) : 

sin^(/xc) = ^^^^ 



The value of at which effects of B{p) become important depends on the value of the quantization 
ambiguity parameter j which following theoretical considerations of Refs. 0,E3|> been set equal to 
half. 
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FIG. 9: A comparison of the quantum dynamics represented by the expectation values and dis- 
persions of \v\^ (points with error bars), the effective dynamics (dotted curve) and the classical 
dynamics (dashed curve) is presented near the bounce point. The effective dynamics provides a 
good approximation to the underlying quantum theory. The values of parameters were = 10000 
and e = 2. 



and using Eq. ()B3|) . It turns out to be 

with p = p'^/2\p\^. Quantum geometry effects thus lead to modification of the Friedmann 
equation at the scales when p becomes comparable to pcrit- The Hubble parameter vanishes 
when energy density becomes equal to the critical value pcrit and the universe bounces from 
the expanding (contracting) branch to the contracting (expanding) branch. Further, p/ p^it 
term becomes negligible for p <^ pcrit and Eq. ()B6|) reduces to the standard Friedmann 
equation; there are no departures from classical general relativity in the low curvature regime. 
Detailed phenomenological investigations of the effective theory obtained from the improved 
quantum constraint confirm that unnatural effects of the old quantization are cured in the 
improved quantization [32]. 

To compare effective dynamics with the exact quantum evolution in Sec. IV Bl it is useful 
to combine ()B3|) and ()B4|) to obtain: 

^ = Vl2^(l-J^Y\. (B7) 

d(p V Pcrit / 

Trajectories obtained from this equation are plotted in FigjU] and compared with the ex- 
pectation values of \v\^ and the classical dynamics. As can be seen, the effective dynamics 
provides an excellent approximation to the underlying quantum dynamics. 



37 



Remark: The modified field equations ()B6|) and ()B7|) can be interpreted as saying that 
the effective Newton's constant is given by Gefr = G(l — p/pcrit), where G is the low energy 
Newton's constant and pcrit ~ 0.82ppi. Now, the renormalization group analysis based on 
Euclidean quantum gravity [s^ strongly suggests the existence of a non-trivial fixed point at 
which the theory becomes asymptotically free. The behavior of Gefr in LQC is in qualitative 
agreement with that picture. 

2. Subtleties associated with = cosmologies 

To write the space-time metric in the standard FRW form, gah = —'Vat'Vbt + a'^it)q°i^, one 
needs a fiducial, Riemannian 3-metric g^^ of constant curvature k. Clearly, under a rescaling 
o?q°ahi the scale factor a{f) scales as a(t) oT^aif). As is well-known, while this 
rescaling freedom can be eliminated in the /c = ±1 cases by requiring that q°^^ be the unit 
3-sphere or unit 3-hyperboloid metric, in the A; = case there is no natural way to select 
a unique whence the scale factor a(t) by itself does not have a direct physical meaning. 
This feature introduces certain subtleties in the Hamiltonian framework in the = case. 
The purpose of this appendix is to summarize them and correct certain misconceptions that 
have permeated in some of the LQG literature. 

The Hamiltonian framework has more information than that contained in the equations 
of motion. In particular the symplectic structure can be regarded as the imprint left by the 
quantum theory on the classical framework. To define the correct symplectic structure and 
Hamiltonian in the non-compact, homogeneous context, one has to introduce a fiducial cell 
V and restrict all integrations to it. In the k = —\ case, one can fix V to have unit volume 
with respect to the unit hyperboloid metric. In the A; = case, such a simple prescription 
is not available. However, once V is chosen, its presence can be used to make the canonical 
variables, the symplectic structure and the Hamiltonian insensitive to the choice of the 
fiducial metric g"^ [4] . This was the procedure followed in as well as in the present paper. 
For each choice of V we obtain a physical theory. In that theory, one can make well-defined 
predictions. However, since in the = case there is no physical procedure to fix the fiducial 
cell, we now have a new 'gauge freedom' V X^V. Under this rescaling the theory changes. 
Only those quantities which are insensitive to this rescaling can be regarded as physical. 

Recall that in the quantum theory V = \p\'^^'^ is the volume operator associated with the 
fiducial cell V. Thus, |//) (respectively \v)) is an eigenstate of quantum geometry in which 

the volume of V is \pf/^ = (Sn-f/G)'^/^ IpI^^^^i (respectively (8717/6)2/2 (3^373/2^2) w^i)- 
Now, in the quantum theory based on the cell V, |l/pp/2 ceases to be a good approximation 

to the eigenvalues of the inverse volume operator 1/V for p < p^, for a well-defined real 
number p^, (which makes no reference to the size fiducial cell V). Therefore one sometimes 
finds statements in the LQC literature to the effect that the quantum geometry effects 
become important when p < p^. However, physically this would mean that the quantum 
geometry effects become important when the physical volume of V is less that Ip^l^^^. Clearly, 
since the fiducial cell cannot be fixed unambiguously in the k = case, unlike in the k = ±1 
cases, this tempting suggestion has no invariant meaning. Put differently, from a space-time 
perspective, one can begin with a quadruplet {qab, 4>', Qab, 4>) representing the initial data 
at late times for an expanding universe, consider a semi-classical quantum state which is 
sharply peaked about this configuration and ask when, in the backward quantum evolution, 
there are significant departures from the classical trajectory. The answer to this physical 
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question cannot depend on the choice of the fiducial cell used in the construction of quantum 
theory. Therefore, in the k = case, the answer that this happens at p = is not viable. 
Similarly, in the analysis of the quantum bounce, the value of Pcrit (or fcrit) at which the 
bounce occurs has only mathematical meaning within the theory based on a specific V. 
From a space-time perspective, this is not surprising: it is well-known that the value of the 
scale factor at which something specific happens is tied to the choice of fiducial metric and 
therefore has no invariant meaning. 

On the other hand, in the classical theory the matter density p has direct physical meaning 
also in the k = case from both the space-time and the phase space perspectives. In the 
space-time perspective, it is given by p = (0)^/2 and makes no reference to the fiducial 
metric. From the phase space perspective it is given by p = p|/2|p|^ and is insensitive to the 
choice of V because under the rescaling V A^V, \^p^ and \p\ A^|p|. Therefore 

it is physically meaningful to ask if there is a critical value of density pcrit at which the 
classical theory becomes inadequate. Numerical simulations show that the LQC dynamics 
used in the literature yields pcrit = ^bcritl = k'/{p^), where k and k' are constants that do 
not depend on V and Pcrit refers to the value of p of at the bounce point. The physical 
meaning of \p\ and p^ = (p depends on the choice of V. Since p itself is gauge invariant, 
i.e., is independent of the choice of the fiducial metric or the fiducial cell V, there is a 
conceptual mismatch.^ By contrast, the quantum dynamics presented in this paper leads to 
p^it ~ 0.82ppi for any state which is semi-classical at late times; in particular, pcrit makes 
no reference to the fiducial cell V. 

To summarize, whether one works with geometrodynamics or path integrals or LQG, 
the introduction of a fiducial cell V is essential in quantization of spatially non-compact 
homogeneous models. V can be used to remove the reference to the fiducial metric from 
the theory. However, its presence introduces a new 'gauge' freedom and one has to exercise 
due care in the physical interpretation of the theory. If we have two theories, the first based 
on a fiducial cell Vi and the second on V2 = A'^Vi, the state |p) (or \v)) of the first theory is 
physically the same as the state |A^p) (respectively \\^v)) of the second theory. Using this 
equivalence, one has to ensure that the final physical results are independent of the choice 
of the initial fiducial cell. 

We will conclude with the discussion of a subtlety. Fix a fiducial cell V and consider 
the resulting quantum theory. Then, for v < v^, the inverse volume corrections become 
important and the functional form of p{p^,p) is significantly different from the classical 
relation p = p'^/2p^. Therefore, if a trajectory were to enter this region, the quantum 
corrected p along it does not increase and may never reach the critical value pcrit- Then 
the argument for a quantum bounce would break down. But a detailed analytical and 
numerical examination shows that this could happen only if the assumption p^ ^ h (in 
c=G=l units) is violated. In a state which is peaked at such a small p^, the uncertainty 
in the Dirac observables becomes comparable to the values of the observables themselves, 
whence the state can not be regarded as semi-classical. Such a state falls outside both the 
detailed quantum evolution studied in the main body of the paper and the realm of effective 



Furthermore, even if one works with the theory obtained by just fixing a fiducial cell from the beginning, 
for states with 3> ft (in the c=G=l units) which are favored by semi-classical considerations, pcrit can 
be quite low, violating the expectation that classical general relativity should be valid at low densities. It 
is this second point that was emphasized in and the main body of this paper. The conceptual mismatch 
is a related but separate drawback of that quantum dynamics. 
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equations developed in Appendix IB 11 
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